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The  programs  in  this  report  are  for  use 
within  the  Department  of  the  Navy,  and  they  are 
presented  without  representation  or  warranty  of 
any  kind. 


A  PROCEDURE  FOR  ESTIMATING  AN  OBJECT'S  POSITION 
BASED  ON  TWO  OR  MORE  BEARINGS  WITH  A 
PROGRAM  FOR  A  TI-59  CALCULATOR 

I.   Introduction 

A  procedure  for  estimating  an  object's  position  with 
bearings  taken  on  or  from  two  or  more  stations  is  developed 
in  Section  IV  of  this  report.    In  the  development  of  the  pro- 
cedure, the  following  things  are  assumed:   The  object  and  the 
stations  are  fixed  on  the  surface  of  a  flat  earth  and  the 
position  of  each  station  is  known.  The  error  in  the  bearing 
taken  on  or  from  a  station  is  a  normal  random  variable  with 
a  known  standard  deviation   e   and  a  mean  of  zero  (if  bias 
exists,  it  is  known  and  removed);  and  station  bearing  errors 
are  independent.   The  user  instructions  for  a  TI-59  program 
to  implement  the  procedure  are  given  in  Section  II,  and  the 
program  listing  is  given  in  Section  III. 

As  an  example  to  illustrate  a  use  of  the  program,  suppose 
bearings  are  taken  on  an  object  from  three  stations  (1,  2  and  3) 
as  illustrated  in  Figure  1.   Also,  suppose  that  the  assumptions 
stated  above  are  satisfied  and  that  an  initial  estimate  of  the 
object's  position  is  made  and  that  it  is  relatively  near  the 
object.   This  assumption  is  discussed  in  Section  IV. 


Station  3 


•  \ Initial  Estimate 


Station  1 


Station  2 


FIGURE  1.   Geometry  for  the  Example 


Let  the  measured  bearings  and  bearing  errors  (standard 
deviations)  be: 


=  35 


e   =  4° 
1 


92  =  351° 


e   =  7° 
e2 


=  131 


e   =  5° 
e3 


And  let  the  ranges  and  bearings  of  the  initial  estimate  be: 


r   =  10,000  meters, 


r„  =  15,000  meters, 


r3  =  12,000  meters, 


3-L  =  38° 


62  =  346° 


33  =  127° 


Use  of  the  position  estimation  program  with  this  data 
gives  a  final  position  estimate  (fix)  determined  by: 

x  =  -512  meters 
y  =  -  75  meters 


where   x   is  its  East-West  distance  and  y   is  its  North-South 
distance  from  the  initial  position  estimate.   The  East-West, 
North-South  xy-coordinate  system  with  its  origin  at  the  initial 
estimate  is  shown  in  Figure  2.   So  the  final  position  estimate 
is  512  meters  to  the  West  and  75  meters  to  the  South  of  the 
initial  position  estimate. 

North 

y 


•  Final  Estimate 


Initial  Estimate 


x   East 


FIGURE  2.   The  Location  of  the  Final  Position  Estimate  with 
Respect  to  the  Initial  Position  Estimate. 
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Minimum  area  elliptical  confidence  regions  for  an  object's 
position  can  also  be  found  by  using  the  Tl-59  program.   The 
centers  of  the  regions  are  at  the  fix,  and  their  axes  lie  along 
the   x'   and   y'   axes  of  the  coordinate  system  obtained  by 
rotating  the  East-West,  North-South  xy-coordinate  system  with 
origin  at  the  fix  through  an  angle   Y .   The  angle  Y   is  defined 
so  that  it  is  positive  for  a  rotation  in  the  counterclockwise 
direction. 

With  the  data  from  the  above  example,  the  program  gives 
Y  =  -31°;  so,  the   x'   axis  is  directed  31°  South  of  East.   For 
a  confidence  region  with  minimum  area  and  a  confidence  level  of 
.9000,  the  ellipse  bounding  the  region  has  a  semi-major  axis 
of  2064  meters,  and  a  semi-minor  axis  of  1453  meters.   The  area 
of  the  region  is  9.4  3  square  killometers  or  2.75  square  nautical 
miles.   The  region  is  shown  in  Figure  3. 


x  East 


FIGURE  3.   A  .9000  Confidence  Region  for  an  Object's  Position 


In  the  example  discussed  above,  the  position  of  the 
initial  estimate  is  an  input  to  the  program.   If  this  is 
not  desirable,  the  program  can  be  used  to  determine  a  po- 
sition for  the  initial  estimate.   The  position  is  the  in- 
tersection of  the  two  bearing  lines  corresponding  to  the 
first  two  bearings  entered  in  the  program.   Both  options 
are  illustrated  in  Section  II. 

Since,  in  general,  the  smaller  the  bearing  errors, 
the  more  likely  that  the  initial  estimate  will  be  relatively 
near  the  object;  small  bearing  errors  can  be  considered  to 
be  a  condition  on  the  use  of  the  procedure. 

Note,  if  the  length  of  the  base  line  joining  the 
first  two  stations  is  small  enough  and  their  bearing  errors 
are  large  enough,  observed  bearing  lines  from  the  two 
stations  may  not  intersect.   If  they  do  not  intersect,  the 
initial  estimate  determined  by  the  program  will  be  at  the 
intersection  of  the  reciprocal  bearing  lines,  and  a  gross 
error  can  result. 


II .   User  Instructions 

The  TI-59  program  to  which  the  user  instructions  in  this 
section  apply  can  be  used  to  calculate  the  quantities  described 
in  Section  I. 

The  program  requires  the  following  inputs: 

1.  the  observed  bearing  from  or  on  an  object  for  two  or  more 
stations; 

2.  station  positions  relative  to  a  reference  position;  and 

3.  the  bearing  error  (standard  deviation)  for  each  observed 
bearing. 

Station  positions  can  be  specified  in  either  of  two  ways. 
In  the  first  way,  Mode  A,  each  station's  position  is  specified 
in  terms  of  its  bearing   a   and  its  range   p   from  a  reference 
position.   In  the  second  way,  Mode  B,  each  station's  position 
is  specified  in  terms  of  its  East-West  distance   x   (plus  for 
East)  and  its  North-South  distance  (plus  for  North)  from  a 
reference  position.   The  reference  position  can  be  any  convenient 
location.    For  example,  if  it  were  at  a  station,  then  for 
that  station   a  =  0   and   p  =  0   or   x  =  0   and   y  =  0. 

The  program  also  requires  an  initial  estimate  of  the  object's 
position.   The  user  has  two  options: 

1.  Let  the  program  provide  an  estimate,  or 

2.  provide  one  with  the  input  data. 

For  Option  l,the  initial  estimate  is  at  the  intersection  of  the 
bearing  lines  determined  by  the  first  two  observed  bearings  entered 
into  the  program.   For  this  reason,  if  this  option  is  chosen,  the 
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first  and  second  groups  of  data  entered  should  correspond  to 
the  two  stations  estimated  to  have  the  smallest  products   r.e.. 
Although  in  this  option  the  reference  position  cannot  be  at 
the  initial  estimate,  it  can  be  at  one  of  the  stations.   If 
only  two  stations  are  involved,  the  final  estimate  is  at  the 
intersections  of  the  bearing  lines.   (If  the  second  option  of 
either  mode  is  used  with  an  initial  estimate  which  is  not  at 
the  intersection  of  the  two  bearing  lines,  the  coordinates  of 
the  final  estimate  will  differ  from  coordinates  of  the  inter- 
section to  the  degree  of  the  approximations  involved  in  the 
estimation  procedure.) 

Two  ways  of  providing  confidence  (probability)  region  data 
are  available.   In  the  first  way,  Mode  C,  the  confidence 
(probability)   p   is  specified.   In  the  second  way,  Mode  D,  the 

multiplier   k   is  specified  where   ka~,    and   ka~.   are  the 

x'  y 

semi-axes  of  the  bounding  ellipse. 

The  values  of  various  quantities  calculated  by  the  program 
are  either  stored  in  registers  or  appear  in  the  display.  If  a 
PC-100A  printer  is  used,  some  of  these  values  will  be  printed. 
The  location  of  calculated  values  and  the  printing  format  is 
given  after  the  user  instructions.  Those  quantities  which  are 
not  described  in  Section  I  are  described  below  in  the  User  In- 
structions or  in  Section  IV. 

All  angles  required  or  calculated  by  the  program  are 
in  decimal  degrees. 


Step  Instructions  Enter   Press    Display 

■"■»  1.     If  the  calculator  has  been  in  use  and 
flags  have  been  set  or  the  memory 
repartitioned,  turn  the  calculator  off 
and  then  on. 

■*   2.     Read  Side  1  and  Side  2  of  Card  1. 
-*3.    Read  Side  3  of  Card  2. 


MODE  A:   Station  Locations  Specified  in 
Terms  of  Bearing  and  Range  from 
a  Reference  Point. 

4a.    If  the  initial  position  estimate  will 
be  determined  by  the  program,  go  to 
Step  7a.  See  the  note  on  Page  10. 

5a.    Enter  the  initial  estimate's  bearing.        a*      A' 

6a.    Enter  the  initial  estimate's  range.  p*      r/s 

7a.    Enter  the  measured  bearing  on  the  0.      a 

object  from  a  station  or  the  reciprocal 
of  the  measured  bearing  on  a  station 
from  the  object. 

8a.    Enter  the  station's  bearing. 

9a.    Enter  the  station's  range. 

10a.    Enter  the  bearing  error. 

11a.  Repeat  Steps  7a,  8a,  9a  and  10a  for  all 
stations.  The  number  of  repetitions  i 
appears  in  the  display  after  Step  10a. 


ai 

R/S 

Pi 

R/S 

e-i 

R/S 

Step  Instructions  Enter   Press      Display 

•*  MODE  B:   Station  Locations  Specified  in  Terms 
of  East-West  Distance  and  North-South 
Distance  from  a  Reference  Point. 


v  4b.    If  the  initial  position  estimate  will 
be  determined  by  the  program,  go  to 
Step  7b.   See  the  note  on  Page  10. 

5b.    Enter  the  initial  estimate's  East-West      x*       B' 
distance. 

6b.    Enter  the  initial  estimate's  North-         y*       R/S 
South  distance. 

^  7b.    Enter  the  measured  bearing  on  the  8^        B 

object  from  a  station  or  the 
reciprocal  of  the  measured  bearing 
on  a  station  from  the  object. 

*8b.    Enter  the  station's  East-West 
distance. 

^9b.    Enter  the  station's  North- South 
distance. 

*10b.    Enter  the  bearing  error. 

*llb.    Repeat  Steps  7b,  8b,  9b  and  10b  for  all 
stations.   The  number  of  repetitions  i 
appears  in  the  display  after  Step  10b. 


BOTH  MODES 

*12.     Calculate  the  East-West  distance,  R/S 

the  North-South  distance,  the 

bearing  and  the  range  of  the 

position  estimate  relative  to  the 

reference  position.   Also  calculate 

the  rotation  angle  y ,  and  the 

standard  deviations  C^,      and  O^,. 

x'        y  * 

To  include  additional  bearing  measure- 
ments after  this  calculation,  go  to 
Step  18. 


xi 

R/S 

y± 

R/S 

ei 

R/S 

Step Instructions Enter Press Display 


13.    For  confidence  (probability)  region 

calculations,  go  to  Step  14  if  the 

confidence  (probability)  for  the 

region  is  specified.   If  k   is 

specified  where   ka^.   and  kaA1 

x'         y  • 

are  the  semi-axes  of  the  bounding  ellipse 

with  the  larger  the  major  axis,  go 

to  Step  16. 

**  14.    Enter  p,  the  confidence  level  p        C        Area 

(probability)  and  calculate  k, 

ka..,  kn   and  the  area  of  the  region. 
x'     y ' 

(The  area  units  correspond  to  the 
distance  units  used.) 

"tfl5.    For  a  different  value  of  p,  go  to 
Step  14. 

16.    Enter  k  and  calculate  the  confidence      k        D        Area 

level  (probability)  p,  kcu.,  ka,,.   and 

x'    y ' 

the  area  of  the  region .   (The  area 

units  correspond  to  the  distance  units 

used.) 


17.  For  a  different  value  of  k,  go  to 
Step  16. 

18.  To  include  an  additional  bearing  measure- 
ment from  either  a  new  or  old  station, 
go  to  Step  7a  if  using  Mode  A  or 

Step  7b  if  using  Mode  B. 


NOTE:   If  a  data  entry  error  occurs  in  either  mode, 
press  RST  and  then  use  the  following  pro- 
cedure:  For  Option  1,  return  to  Step  7 
and  repeat  all  data  entries.   For  Option  2, 
return  to  Step  5  and  repeat  all  data  entries, 

Also,  if  a  position  estimate  is  to  be  de- 
termined for  a  new  object  position  or  if 
a  new  mode  is  to  be  used,  follow  this 
instruction. 
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NOTES : 

a)   The  program  printing  format  is  given  below: 

For  the  initial  data,  i  =  1,  2,  ...  ,  n  with  one  space 

between  groups: 


Mode  A 


a" 


9i 


a. 


initial  estimate 
if  provided 


Mode  B 


e. 

1 


e. 

l 


e . 

i 


The  format  for  the  calculated  position  data  is 

x 

y 

p 

Y 


y' 


li 


For  the  confidence  (probability)  region  portion  of  the  program 
the  format  is  after  pressing  either  C   or  D: 

P 
k 

ko~ ,  semi-axis 

x1 

ko ~  ,  semi-axis 

y' 


Area 


b)   The  following  data  is  stored  in  the  indicated  registers: 

Data  Registers 

x*  R38 

y*  R39 


Y 


R2  9 


x  R30 

y  R31 

a  R32 

P  R33 


a  -  R16 

x ' 


a  a  R17 

y' 


p 

R14 

k 

R15 

ka~  . 
x' 

R18 

y' 

12 

R19 

Four  data  tapes  for  a  sample  problem  are  given  below. 
Distance  units  have  not  been  specified,  but  they  could  be 
meters  for  example.   Angles  are  in  degrees.   Option  1  (initial 
estimate  not  provided)  for  Mode  A  and  Mode  B  is  indicated  by 
A  and  by  B  and  Option  2  (initial  estimate  provided)  is  indicated 
by  A ■  and  B ' . 

For  each  mode  and  each  option,  the  input  data  are 
indicated.    The  data  determine  the  relative  locations  of  three 
stations  as  well  as  the  observed  bearing  of  an  object  from 
each  station. 

For  A  and  B,  the  reference  location  is  at  Station  1 
and  the  initial  position  estimate  (determined  by  the  program) 
is  at  the  intersection  of  the  bearing  lines  for  Station  1  and 
Station  2. 

The  intersection  has  coordinates   x*  =  906.4853528 
and   y*  =  17296.77092   with  respect  to  Station  1. 

For  A*  and  B',  both  the  initial  estimate  and  the 
reference  location  are  at  the  intersection  of  the  bearing  lines f 
so   a*  =  0   and   p*  =  0   and   x*  =  0   and   y*  =  0. 

The  data  for  A,  B,  A*  and  B'  are  all  equivalent, 
and  each  solution  gives  the  same  data  for  a  confidence  (prob- 
ability) region  calculation.   A  tape  with  confidence  (prob- 
ability) region  results  for  both  Mode  C  and  Mode  D  which 
correspond  to  A,  B,  A'  and  B'  is  given  with  the  first  four 
data  tapes. 
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u. 

4 


10000. 

303. 
14000, 


°1 

al 

Pi 
e. 


°2 

a2 

P2 


e3 

0U 


.  '  .   ''-:  '   .••  :■ ;  '-H  :  : 

X 

6462. 71223 

Y 

. 995471725 

a 

6  4  7  2 .  7"  0157 

P 

, 325392245 

Y 

C;  ~   "'  i  £  "'  y  cr  s." 

a 

0. 
0. 


320. 50808 


4. 


a* 
P* 


ei 

ai 
pl 


D2 
a2 

20000.      P2 
3.       e2 


-J  LI  3. 


1  •-  q      =• -i  .•■  -?  -j  tr  c 


x' 


!97836     p 


3':!id«  8974567  x 

834. 0536335  v 

20 1. 7534019  a 

893. 03931 1 i  P 


i 


325392259 


i  O  ■'  .    O  C'  C'  •-<  (    J  i' 

1 2  3  3 .  0  8 0776 


Y 
a  „ 


y' 


14 


B 


B 


C    or   D 


0. 
0. 
4. 


-99Ste 


. 295348 


82 

x2 


523. 3595624     y2 


303.      9  3 
7624. 94649    x 
11741.38795     y3 


573. 5878927  x 

16462. 71223  y 

1.995471723  a 

16472.70157  p 


i   ■   •_■  i-  -J  •_•  7  £.  i-  "t  >J       I  - 

787.  3663756    <?-, 
-___    „x' 


1233.  080' 


y' 


0.       x* 

U.  y* 


Hli-"  4*53528      Xl 

1  .-'296.  .','092      y1 

4.        e. 


33.  e2 
-10892.  7807  x2 
16773.41136      y. 


303.  63 

6718.  461 137  X3 

•5555.  382969  y 

r. "  3 

'-•  e3 

■332.  8974589  x 

•834.  0586906  y 

201. 7584019  a 

898.0393184  p 


7. 325392239 
787. 3663755 
1233. 080777 


0.  9 

P 

*»  .-;  c."  -•*■•  •  "  ,-. .— .  .- 

i  4  D  y  h  o  U  £  b 

k 

89.  t-  1493 

ka~, 

46, 149454 

!    X 

046364.  97 

Area 

b  b  4 


i  •- '  i  *f  i   ■_',_.  ._!  1 

ka 

2466. 161553 

ka|: 

1  22005  i  7.  6 

,  y 

Area 

x' 
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The  following  data  tape  illustrates  the  effects  of 
using  only  the  bearings  for  the  first  two  stations.   The  tape 
is  for  Mode  A,  Option  1.   The  values  stored  in  Registers  38 
and  39  (the   x  and  y  coordinates  of  the  intersection  of  the 
bearing  lines  from  Station  1  and  Station  2  with  reference, 
in  this  case,  to  Station  1)  are  also  listed  on  the  tape  (as 
well  as  given  above).   And,  as  can  be  seen,  the  initial 
estimate  and  final  estimate  correspond. 

The  data  tape  also  illustrates  the  use  of  additional 
bearing  data  to  revise  a  position  estimate.   The  data  for 
Station  3  printed  after  the  first  confidence  (probability) 
region  calculation  results  was  entered  by  again  repeating 
Step  11a,  and  the  remaining  results  were  obtained  by  next 
repeating  Step  12  and  then  Steps  14  and  16.   Note,  these 
results  are  the  same  as  the  corresponding  results  for  Mode  A 
on  page  14. 
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To  obtain  the  results  given  in  Section  I,  use  A' 
and  take  the  reference  position  at  the  initial  estimate 
(a*  =  0,  p*  =  0).   Then   (*-,_  =  218°,   a2  =  166°   and 
a_  =  307°.   The  data  tape  for  the  calculation  is  given  below 
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III.   Program  Listing 

Before  entering  the  program,  press  2nd  and  then  CP  or 
turn  the  calculator  off  and  then  on.   Next  enter  5  in  the  display, 
press  2nd  and  then  Op  17.   This  repartitions  the  calculator's 
memory  so  that  the  complete  program  can  be  entered. 

Before  recording  the  program,  enter  6  in  the  display, 
press  2nd  and  then  Op  17.   This  returns  the  calculator's  memory 
to  the  normal  partition  (479.59).   Returning  the  calculator  to 
the  normal  partition  allows  the  two  program  cards  to  be  read 
in  the  normal  partition  without  forcing.   When  the  program 
is  used,  it  repartitions  the  calculator  so  that  Bank  3 
registers  are  program  registers. 
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IV.   A  Development  for  the  Procedure 

In  the  development  for  the  estimation  procedure  given  here, 
all  angles  are  in  radians  and  the  assumptions  stated  in  Section  I 
apply. 

Figure  4  shows  three  bearing  lines  from  the  ith  of   n   stations. 
One  is  the  observed  bearing  line  of  an  object.   One  of  length  K.      goes 


Initial 
Estimate 


forth 


x 


Observed 
Bearing 
Line 


to  the  origin  of  an  xy-coordinate 
system  located  at  the  object's  un- 
known position.   And  one  of  length 
r.   goes  to  an  initial  estimate 
with  known  position  but  unknown 
coordinates   (x,y) .   Note,  esti- 
mates for   -x   and   -y   estimate 
the  object's  position.   To  find 
estimates   -x   and   -y,  consider 

the  arc  coordinates   u.  =  >t .  (9  . -<J> .  ) 

1     11  ri 

of  the  observed  bearing  line  and 

v.  =  1 . ( $ . -6 . )   of  the  bearing  line 
l     11  ri  r 

to  the  point   (x,y) .   They  are  de- 
fined by  the  three  bearing  lines 
and  the  circle  of  radius 
goes  through  the  object's  position 
and  which  is  centered  on  the  station 
as  shown  in  Figure  4. 


K  .      which 

i 


Station  i 


FIGURE  4.   Problem  Geometry. 
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By   defining     w.    =   ^ (8 . -6 . )       (all    angles    in   radians), 

u.    =  v.    +  w..      Note,       6.-3-       is   known,    but      3  •-<{>■       is   not. 
ill  'ii  ii 

However,    v.       can  be   expressed   in   terms   of      x      and      y,    and,    to    first 
order,    v.    =   x   cos    8.    -   y   sin    e .  ;    so,    if     tan    (3--4>j)    -    (B  ■-<()■)       for 

1  1  1  X     X  X     X 

i  =  l,2,...,n,  that  is,  if   (x,y)   is  relatively  near  the  object's 
position,  u.  =■  4.(6.-$.)  +  x  cos  8.  -  y  sin  6.   for   i  =  l,2,...,n. 

In  this  development,  the  initial  estimate  is  assumed  to  be 
relatively  close  enough  to  the  object's  position  so  that  the  above 
approximation  for   u.   can  be  used  and  so  that  H..      can  be  taken 
equal  to   r..   With  this  assumption,   u.   is  a  function  of:   9., 
the  observed  value  of  a  random  quantity;  the  known  parameters   r. 
and   3-;  and  the  unknown  parameters   x   and   y. 

If  a  distribution  for  the   0 .   can  be  specified,  then  a 

distribution  for  the   U.   can  be  determined  and  standard  estimates 

l 

x   and   y   for   x   and   y   can  be  considered.   Maximum  likelihood 

estimates  are  discussed  in  this  section.     Each   ©.   is  taken  to 

2 
be  a  normal  random  variable  with  mean   <j>  .   and  variance   e..   And 

i  l 

the   n   random  variables   0.,  i  =  1,2,.. . ,n   (one  for  each  station) 
are  taken  to  be  independent. 

The  likelihood  for  a  sample   6.,e_,  ...  ,  6    is  then 

r      12         n 

L(ei'e2 v  =  l  (-^tr)  exp  t"  1 1  (w2/ei] 

and  the  likelihood  for  a  corresponding  sample   u,  ,  u„ , . . . ,  u   is 
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(Ul,u2,...,un)    =   ^  (_-i_)exp  [-  \  |   u^] 


where      a.    =   r.e.       (with      e.       in   radians)    since      u.    =   r.(G.-d).). 
111  1  111    vi 

By  definition,  the  maximum  likelihood  estimates  of   x 

and   y   are  the  estimates   x   and   y  which  make   L (u, ,u„ , . . . ,u  ) 

1   2      n 

a  maximum.   In  this  case,  making   L (u, ,u» , . . . ,u  )   a  maximum  is 

L      2  n 

equivalent  to  making   £, (u./a.)   a  minimum.   So,  to  find   x 
and   y,  solve  the  following  two  equations  for   x   and   y: 


9(£n  L)        a,     3(&n  L) 

=  U    and    ~ =  0  . 

3x  3y 

The  solutions  are   x  =  x   and   y  =  y,    and   x   and   y   are  the 

maximum  likelihood  estimates.   With   w.  =  r. (0.-3.)   and 

l     ill 

the  conditions  assumed  above  these  two  equations  are  linear 
equations  in   x   and   y.   And, 


n  2 

]  (cos  3. )/a   =  0 


)  [w.  +  x  cos  3-  -  y  sin  3- 
k        i  i   J  i 


and 

n 


)  [w.  +  x  cos  3-  -  y  sin  3-]  (sin  3-)/cr.  =  0 
l         -l  ^i         J  1        i  '     i 


And,  in  terms  of  the  following  quantities: 

2       2  2 

A  =  Z(cos   3-)/a.  ,      B  =  Z(sin  3.  cos  3 ■ ) /a .  , 


C  =  Z(sin2  3-)/a2  ,      D  =  Z(w.  cos  3-)/a2  , 
li  i       li 


E  =  E(w.  sin  3- )/a2  , 
i      11 
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the  equations  are : 


So  the  solutions  are: 


Ax  -  By  =  -D 
Bx  -  Cy  =  -E  . 


x  =  (BE  -  CD) /(AC  -  B2) 
y  =  (AE  -  BD)/(AC  -  B2) 


A  confidence  region  can  be  constructed  about  an  estimated 
position.   In  order  to  indicate  how  this  is  done,  a  probability 
region  about  the  true  position  will  be  considered  first. 

Note,  x   and   y   are  values  of  random  variables.   If  a 
new  set  of  bearings   9, , 0~  ,  .  .  .  ,  0    is  observed  (for  a  fixed  initial 
estimate  and  object) ,  in  general,  a  new  pair  of  values   x   and   y 
will  be  obtained. 

If   X   and   Y   represent  these  random  variables,  then 

-i     n       2 

X  =  i-^—  I    (W./af)  (B  sin  6.  -  C  cos  3.) 

(AC-B^)  1    *   1  1  1 

Y  =  ^—  I    (W./a.)  (A  sin  3.  -  B  cos  3.) 

(AC-B2)  1    x      1         1  1 


with   W.  =  r. (0.-3-).   (W.   is  the  random  distance  intercepted 
l     11   l       l  e 

along  the  ith  arc  between  the  bearing  lines  defined  by   0.   and 
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Note,  X   and   Y  have  a  bivariate  normal  distribution, 

since  they  are  a  linear  combination  of  the   n   normal  random 

variables   W..  ,W_,  .  .  .  ,W  ,  or  equivalently  of  the   n   normal 

random  variables  ©  ,© ~,  .  .  .  ,©  .   Also   E  (W.  )  =  r .  (d> .  -3  .  )  . 

12'    '  n  l      i  Yi   i 

If   3-  =  <J> .   for   i  =  1,2,..., n,  that  is,  if  the  initial 
estimate  of  the  object's  position  is  at  the  object's  position, 
E(W,)  =  0   for   i  =  1,2, ...,n.   And,  therefore,  E(X)  =0   and 
E(Y)  =0.   So,  in  this  case,  the  "location"  of  the  bivariate 
normal  distribution  of  a  point    (X,  Y) ,  the  random 
coordinates  of  the  object's  estimated  position,  is  the  same  as 
that  for  the  point   (-X,   -Y)   and  both  are  centered  on  the 
object's  position.   However,  the  "location"  of  the  distribution 
of   (-X,-Y)   is  independent  of  the  location  of  the  initial 
estimate  when  the  coordinates   (-X,-Y)   refer  to  a  coordinate 
system  with  origin  at  the  initial  estimate.   This  fact  simplifies 
the  establishment  of  a  confidence  region  about  the  location  of 
an  estimated  position. 

A  region  of  minimum  area  for  a  given  probability  of  con- 
tainment of  an  estimated  position  can  be  determined.   The  region 
is  bounded  by  an  ellipse  which  is  centered  on  the  object's  position 
and  whose  axes  lie  along  the  axes  of  an  x'y ' -coordinate  system 
obtained  by  rotating  the  xy-coordinate  system  centered  on  the 

object's  position  through  an  angle  y    .   In  this  system,  o ^  =  0, 

xy 
that  is,  X'   and  Y'   are  independent  normal  random  variables. 
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The  two  coordinate  systems  are  illustrated  in  Figure  5 . 
The  coordinates  of  a  point  in  the  two  systems  are  related  by 
y'  y  x1  =   x  cos  y  +  Y  sin  y 

y1  =  -x  sin  y  +  y  cos  y 

These  relations,  along  with 

o^.^      =  0,  imply: 
x'y" 


FIGURE  5.   Rotation  Geometry. 


x 


2     2  2     2 

a^  cos   y  +  2a    cos  y  sin  y  +  a^  sin   y  , 

x  xy  y 


and 


2  2     2 

a   sin   y  -  2a    cos  y  sin  y  +  a   cos   y 

x  xy  y 


2a 


tan  2y 


xy_ 


2     2 

a   -  a 


x 


where   y,  the  angle  of  rotation  of  the  coordinate  axes,  is 
positive  in  the  counterclockwise  direction. 

With  the  initial  estimate  of  the  object's  position  at 


the  ob 


ject's  position   (3.  =  <j> .  ,  i  =  l,2,...,n),  so   E(W.)  =  0 


and  Var (W. )  =  a. , 
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2 

a  *    =   - 


x    (AC 


.2 


i-T-T  I  (Va?)  <B  sin  6,  -  C  cos  3.)2, 
-BZ)Z  1      1  X  1 


5-^-  J  (1/a2) (A  sin  6.  -  B  cos  3.) 

-B^)z  1     1         1  x 


y    (AC- 


and 


a 


/N  >N 


1     n      2 
— 5-^5-  £  (1/a .  )  (B  sin  3  .  -  C  cos  3  .  )  (A  sin   6  .  -  B  cos  6  . 

-B  )   1  1  -1  1  1 


xy    (AC- 
Using  the  definition  for  A,  B   and  C,  the  above  become 


2    _C 

x    (AC-BZ) 


a2  =  -± 


(AC-B2) 


and 

B 


xy    (AC-B2) 

So,  tan  2y  =  2B/(C-A)   for   3.  =  4>  •  ,  i  =  1,2, ...,n. 

With  the  object's  position  known  and,  hence,   <J> .   known 

2    2 
for   1  =  1,2,. ..,n,  the  above  equations  for  a    .    o    ,    a  and   y 

A  A  A  A 

x   y   xy 
can  be  used,  since  the  initial  estimate  of  the  object's  position 

can  be  taken  as  the  object's  position. 

With  values  for   aA,  aA,  aAA   and   Y,  values  for   aA 

x   y   xy  x' 

and  a*        can  be  found  by  using  the  equations  in  the  middle  of 

r 

Page  30.  And  then,  the  probability  that  an  estimated 

position  will  be  within  an  ellipse  of  semiaxes   kaA   and  kaA 

x'         y' 
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which  is  centered  on  the  object's  position  can  be  found.   It  is 

2 

i  -  exp(-k  /2)  .   (This  result  follows  from  integrating  the  bi- 

variate  normal  density  over  the  ellipse.)   And  the  area  of  the 

ellipse   is      irk    aA    o~ 

x1    y ' 

Given  estimates   x   and   y   found  by  using  the  relations 


on  Page  28,  an  ellipse  with  semi-axes   koA    and  k<j A   centered 

x'  y' 

on  the  point  with  coordinates   (-x,-y)   in  a  coordinate  system 

with  origin  at  the  initial  estimate  and  oriented  as  indicated 

2 
by   y   is  a   1  -  exp(-k  /2)   confidence  region.   This  follows 

from  the  bivariate  normal  distribution  of   -X  and   -Y  which 

in  this  system  is  centered  on  the  object's  position.   The  ellipse 

2    2 

is  defined  if  Q^,    o^      and  a^  are  known  (the  covariance  matrix 

x   y       xy 
is  known) .   And  to  the  degree  of  the  approximations  involved, 

this  can  be  assumed  to  be  the  case.   In  particular,  by  assuming 

the  initial  estimate  of  the  object's  position  is  at  the  object's 

position,  which  is  consistent  with  assuming   ($.-<J).)   is  small, 

2     2 

values  for   cr^  ,  a/  ,  a^   and   Y   can  be  obtained  by  using 

x    y    xy 

the  relations  on  Page  31.   These  values  can  then  be  used  to 

2         2 
determine  a^        and  o^        by  using  the  relations  on  Page  30. 

x'         y' 
And,  then,  with  a  value  for   k,  a  confidence  region  can  be 

constructed.   To  the  degree  of  the  approximations  involved,  the 
shape  of  the  confidence  region  is  independent  of  both  the  object's 
position  and  of  the  initial  estimate  of  the  object's  position. 
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For  the  case  where  bearings  are  taken  from  the  object 
on  two  or  more  stations,  e.   is  the  reciprocal  of  the 
bearing  taken  from  the  object. 

A  discussion  for  this  and  for  other  bearings  only 
position  estimation  procedures  for  situations  similar  to  the 
one  considered  here  is  given  in  Reference  1  listed  below. 
Reference  2  gives  an  equivalent  bearings  only  procedure.   It 
also  gives  a  range  only  procedure,  a  range  and  bearing  procedure 
and  HP-9830A  programs  with  which  to  implement  the  procedures. 
Using  the  fix  determined  by  two  lines  of  bearing  as  the  initial 
estimate  was  suggested  by  this  reference. 


The  equations  used  in  the  program  to  determine   (x*,  y*) , 
the  coordinates  of  the  fix,  are: 

x*  sin  (e^e^  =  [  p1  sin  (a^-8^)  ]  sin  62 

-  [ P2  sin  (a2-62) ]  sin  8^ 
y*  sin  (02-9^  =  [  p1  sin  (o^-e^  ]  cos  62 

-  [ p2  sin  (a2-62) ]  cos  Q1 
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